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BUCKLING IK SHEAR OF CONTINUOUS FLAT PLATES 
By Bernard Budiansky, Robert W. Connor, and Manuel Stein 


SUMMARY 


As basic information for the design of thin-veb spars, the 
theoretical shear buckling stress is presented for an infinitely long, 
clamped plate divided into square panels by rigid intermediate supports 
In addition, results are given for the shear buckling stress of a plate 
of infinite length and width having intermediate rigid supports that 
form an array of square panels. The results indicate the fallacy of 
the usually made assumption that each panel buckles in shear as if it 
were simply supported along the intermediate supports. 


INTRODUCTION 


The design of a thin-web beam requires a knowledge of the shear 
buckling stress of the web. In many cases the beam has fairly heavy 
flanges, and many equally spaced, intermediate stiffeners of high 
flexural rigidity and low torsional rigidity. The shear buckling of 
the web of such a beam may be analyzed by considering the web to be 
a long clamped plate, divided into equal panels by intermediate 
supports which completely restrain the plate from deflecting but 
offer no torsional restraint. (See fig. 1.) 

The assumption is implicit in the work of many writers that each 
panel of the plate shown in figure 1 would buckle in shear as if it 
were simply supported along the intermediate supports. (See 
references 1, 2, and 3 -) Although the assumption is valid for the 
compressive buckling of continuous plates having equal bays, it does 
not hold true for shear buckling. Adjacent panels buckling in shear 
restrain each other, so that the plate bending moments do not vanish 
along the intermediate supports, as they must for simple support. 

This continuity effect is probably a maximum for nearly square panels; 
from physical considerations it is evident that the continuity effect 
disappears for very small and very large values of the ratio b/a. 


The present paper gives the theoretical shear buckling stress for 
an infinitely long, clamped plate divided into square bays by inter- 
mediate rigid supports. In addition, as further information on the 
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effect of plate continuity on shear "buckling stresses, results sire 
given for the shear "buckling stress of a plate of infinite length and 
width having intermediate rigid supports which form an array of square 
panels. (See fig. 2.) The theoretical, analyses of "both problems, 
performed by the Lagrangian multiplier method, are presented in 
appendixes. 


SYMBOLS 


a 

b 

t 

E 

D 

T 


length of panel 

width of panel 

plate thickness 

Young* s modulus of elasticity 

Poisson’s ratio 


plate stiffness in bending 


Et 3 

,12(1 -n 2 V 


critical shear stress 

critical-shear— stress coefficient in the 
formula t = k e 



x plate coordinate parallel to length 

y plate coordinate parallel to width 

„ w deflection normal to plane of -the plate 

a mn » ^ran * c mn * ^mn Fourier coefficients 

Lagrangian multipliers 


V internal bending energy 

T external work of applied stress 

m,n,i,j,p,q. integers 

8 m Kronecker delta: 1 if m = nj 0 if m ^ n 
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RESULTS .AND DISCUSSION 

Infinitely Long Clamped Plate with. Square Bays 

* 

The critical shear stress of the plate shown. in figure 1 is given 
"by the formula 


r = k s 


The theoretical analysis given in appendix A yields the result 


k = 13.1k 

s 


for the case of square "bays . This result, which represents the average 
of upper— limit and lower— limit solutions, is within 0.6 percent of the 
true "buckling stress coefficient. 

It Is of interest to compare the present result for a continuous 
plate with the shear "buckling stress coefficients of square plates with 

(a) two opposite edges clamped and the other two simply supported, and 

(b) nil four edges clamped. This comparison Is shown in the following 
table : 



fUpper limit (reference 1). 

“Reference k. 

°Maximum error 0.6 percent (reference 5* appendix A). 


Thus, for square "bays, changing the "boundary conditions a long the 

transverse edges from simple support to continuity provides about 

25 percent of the increase in shear buckling stress that would be 
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obtained, by changing the boundary conditions from simple support to 
clamped edges . 


Doubly Infinite Array of Square Panels 

The buckling in shear of the doubly infinite array of square 
panels continuous over rigid supports shown in figure 2 is analyzed 
in appendix B . The result for the critical— 3 hear-s tress coefficient, 
obtained to four-figure accuracy by means of upper-limit and lower- 
limit solutions, is 


k 0 = 11.10 


Comparison of this result with those for square plates having other 
boundary conditions is shown in the following table: 


Boundary conditions 

k s 

All edges simply supported 

a 9-35 

All edges continuous 

11.10 

All edges clamped 

14 .71 


Reference 6. 


Thus, continuity of the square panels at all edges provides 
over 30 percent of the Increase in shear buckling stress that would 
be obtained by clamping all edges of a simply supported square plate . 


CONCLUDING REMARKS 


The theoretical results presented indicate the fallacy of the 
usually made assumption that continuous plates having equal, finite 
bays buckle in shear as if each bay were simply supported at the 
intermediate supports. The increase in buckling stress due to- 
continuity at the supported edges of square panels is of the order 
of 25 percent of the increase that would-be provided by clamping tiie 
edges. This continuity effect is probably a maximum for nearly square 
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panels; from physical considerations it is evident that the continuity 
effect disappears for very small and very large values of the width- 
length ratio b/a. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., January 20, 19^8 
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APPENDIX A 


INFINITELY LONG CLAMPED PLATE WITH SQUARE BAXS 


The shear hackling of the plate shown in figure 1 will he analyzed 
hy the Lagrangian multiplier method (references 7 and- 5) to obtain both 
upper and lower limits to the true buckling stress. Several possible 
types of buckling configurations will be investigated; the correct 
configuration is that which corresponds to the lowest buckling load. 

Buckling configurations .— It is intuitively evident- that, for 
square bays, each bay will have the same buckling configuration; 
however, it is not evident whether all bays will buckle in the Bame 
direction, or whether the directions will alternate . It is probable 
that the deflection of each square bay of the continuous plate iB 
symmetrical about the bay midpoint since simply supported square plates 
(reference 6) and clamped square plates (reference 8) buckle in shear 
with symmetrical patterns. However, for conpleteness, the possibility 
of an antisymmetrical deflection in each bay will be investigated. 

Boundary and continuity conditions, and deflection functions .— 

From the many possible types of Fourier series discussed in appendix B 
of reference 5j suitable series will be chosen to represent the various 
possible types of buckling patterns. Schematic representations of the 
three buckling configurations considered (where the notations S and A 
refer to buckling patterns that are symmetrical and antisymmetrical about 
the bay midpoint, respectively, and the signs preceding these notations 
indicate direction of buckling), with the series chosen for these 
configurations, and the conditions of continuity which they satisfy, 
are as follows! 



2 xmtx 


m=l n=l 


00 00 
bi= 0 n=0 


2 mmx 


~(o,y) = ^(^y) 
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The conditions of continuity are satisfied "by these deflection 
functions term "by term; but the boundary conditions at the clamped 
edges , 

v(x,0) * w(x,b) » 0 ( A2) 

3 0 (A3) 

and the condition of zero deflection at the intermediate supports, 

v(0,y) = w(b,y) = 0 (A4) 


vill be introduced by means of Lagrangian multipliers. 

The remainder of the derivation will be performed for the first 
of the buckling configurations (Ala) (which, as will be shown, is 
actually the governing one); the details of analysis for the other two 
buckling patterns are analogous to those presented. 

Energy expressions .— The internal energy of the plate V and the 
external work of the stresses T are given by the expressions 




- 2(1 - n) 


Sx 2 5y 2 




5w dw 
3x 3y 


dx dy 


Substituting the chosen expansion for w (equation (Ala)) into 
these energy integrals gives 
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v . ^ Y Y (m 2 + n 2 ) 2 L 2 (l - v, - Son) - ^ 2 (l ♦ 8^ + S 0n ) 
13 ^=6 n=0 L J 

•0 00 

T = 2rtir 2 mna^d^ - 

m=l n=l 

Then 



where 

A^ = 2(m 2 + n 2 )(l + 5^ + 5^) 

Note that V — T is independent of clgg, since Aqq = 0. 

Constraining relationships In order to satisfy the boundary 
conditions of zero deflection (equations (A2) and (Ah)), it is necessary 
to impose the following constraint relationships: 

^ djjjj =0 (j = 0jlj2,..») (A6a) 

111=0 

00 

d in = 0 (i = 0,1,2,...) (A6h) 

n=0 
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Similarly, in order to satisfy the condition of zero elope along 
the clamped edges (equation (A3)), it must he true that 


^ na in " 0 (i = 1,2,3,...) (A7) 

n=l 


The term d^Q, which does not appear in the energy expression V — T 

(equation (A 5 )), may he eliminated from these constraining relationships 
hy subtracting the first of equations (A6b), the equation for i = 0 , 
from the first- of equations (A6a), the equation for J = 0 . The necessary 
constraining relations finally become 


CO CO 

Z*--Z 

m=l n=l 


d 0n = 0 


J — ® ( <3 = 1,2,3,.*») 


(A8a) 


m=0 


^ d in = ® ( d “ 3 -> 2 , 3 ,«»») 


n=0 


CO 

^ na^^ =0 (i =1,2,3,...) 


(A8b) 


n=l 


Lower— limit solution .— A lower— limit solution requires that V — T 
(equation (A5)) he minimized with respect to all the coefficients a^ 

and 4 ^, whereas the constraining relationships (equations (A8a) 

and (A8b)) are to he satisfied only as far as i = -p and J = q. As 
required hy the Lagrangian multiplier method, the expression to be 
minimized is 
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k s mna mn d ran ( 


- a, 


m=l 


d m0 ~ 'y X i ‘Vj 

j=l m =0 


n=l 



(A9) 


e c h t .„ >r - 


The quantities a, Aj, ji^'and T) i are Lagrangian multipliers. The 

equations for minimizing V — T while satisfying the constraining 
relationships on the a’s and d*s (equations (A 8 a) and (A 8 b) up 
to i = p and J = q) become 


* =0 
3a mn 

(m,n “ 1 ^ 2 , 3 * •••) 

90 =0 
^mn 

(m,n = 0 , 1 , 2 , . ..) 


Equations (A 8 a) and (A 8 b) up to i = p and J = q 


(A 10 ) 


By evaluation. 


as 

^mn 


“ ^mn^han ~ ^mrid^ nT toi 


= 0 


(All) 


as 


(A12) 


= Amn'W ~ kg^ne^ “ x n “ ^ = 0 ( m * n f 0) 
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5J“ = “ a “ hn = 0 (A13) 

= ^011^011 + a " *n = ® (Al^) 


In equations (All) to (Alt) X n does not appear for n > q, nor 
do T) m and for m > p. When both m > p and n > q, the 

multipliers \ n , ri m , and n m vanish from equations (All) and (A12). 
Then one of two conditions is possible: either 

p 2 2 2 

- k g in = 0 
or 

a mn = <i mn = ® 


The first alternative, however, for given values of m > p 
and n > q, will ordinarily lead to a very high value of the buckling 
coefficient k 0 . For the lowest buckling load, therefore, when m > p 

and n > q. 


a mn “ ^mn = ® 

For the remaining a*s and d’s-, solving equations (All) and (A12) 

gives 

“ ^mn^m + c mn( x n + ^hn) ^a) 


d mn = B mn(*n + hn) + nC mn T bii 


(A15b) 
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where 


®mn 




^Tm 




'mx 


x m ? - »m? 


in which 


D mn = k s m 


From equations (AI3) and (AlU), 


d ^>'i( a + u ») <Al6) 

d 0n “ A^(" a + x n) (AJ.7) 


It, is to he emphasized that, in equations (A15a) to (A17), 


and 


= u m = ° for m > p 
X n = 0 for n > q 


Substituting the values of and d^ given by equations (A15a) 

to (A17) in constraining relationships (A8a) and (A8b) up to J = q 
and i = p gives 



oo \ p q 

a ( XI Bmd + 2Z BQn ) + XI - X! Bor ^ n " 0 

,m=l n=l / m=l • n=l 


(Al8a) 


H 

•f 


^)Y_ 3 m)-^ 0 3 * 2l V* + JL J! X)V * 0 (J * l-M'— 1 > 

m=0 m=l jn=l 


(Al8b) 


Hi 


°° q co 

X B ln + ^10 + > B ji^n + 

n=0 n=l n=l 


nC ln » 0 (i = l,2,3,...p) 


(Al8c) 


h'^ r ‘\n + li 1 !£ in\ + '' l '^'* : 1 j t '- 0 (1 = 1 , 2 , 3 ,... P) 

n=l n=»l n*=l 


(Al8d) 


In order for these 2p + q + 1 equations to he compatible, the determinant of the coefficients of 
the Lagrangian Multipliers must vanish. From this determinant*! equation the critical value of the 
buckling coefficient may be found. Those terms in the stability determinant involving infinite 
summations may be evaluated by means of computation aids similar to those described in the appendix 
of reference 8. 

Upper-limit solution .— The theory of the upper-limit solution in the Lagrangian multiplier 
method (reference 7) requires that some a's and d*s arbitrarily be set equal to zero, that 
equation (A5) be minimized vith respect to all the remaining a's and d's, and that all the 
constraining relationships of equations (Afia) and (A8b) be satisfied. 
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As a result of the necessity for satisfying all the constraining 
relationships, a redundancy exists among equations (A8a) (see 
reference 7)j this redundancy can he removed hy discarding the first 
of equations (A8a) . The necessary constraint relationships now become 


X^-° 

m=0 


OO 

X 4 in ■ 0 

n=0 


m 

y~ “in = 0 
n=l 


The arrays of Fourier coefficients to he retained in an upper- 
limit solution are chosen with the intention of combining accuracy with 
ease of solution. The procedure followed in the analysis of the shear 
buckling of . a clanged plate (reference 5, appendix A) suggests that the 
restriction on the existence of the coefficients can be as follows: 


(J - 1,2,3,...) (A 19 a) 


(i - 1,2,3,...) (A19b) 


(i - 1,2,3,...) (A19c) 


d |U| =* 0 (when either m > p or n > q) 


ajjQ = 0 (when a > p) 


When these limits are imposed on the existence of the coefficients, 
the constraint relationships (equations (A19)) take the form 

P 

X = 0 (J = 1,2,3, ...q) " (A 20 a) 

m=0 


<1 

X ^■in “ ® (i “ l,2,3,..«p) (A20b) 

n=0 
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n=l 


( i - 1, 2j 3 > • • »p) 


( A20c ) 


The function to he minimized is 


G = . Iv Ln 2 ^ - SfflD - Son) + ^mn 2 ] - k S mna mn d mn f 

m=0 n=0 L. J 


J=1 m=0 i=l n=0 i=l n=l 


na 


in 


(A21) 


Setting = 0 then gives 

a mn " + c mn( x, n + (for / 0) 

a mn = ( for = °) 

jtnn 

V = B mn( x n + ^m) + n( Wlm ( for m,n ^ 0) 
***> = A^ m 

= A^n 


Substituting these values in the constraint relationships (equations (A20)) 
gives the final stability equations: 



I 


p p p 

+ y ~ B mA + y~ * 0 (J - 1,2,3,... l) (A22a) 

tt=«0 m=»l m=l 


q. <1 q. 

B ir^n + 51 B ln + ^ln “ 0 (i - 1,2,3, ...p) (A22b) 

n°l n=0 n=l 




+ 



n«l 


0 


(i = 1,2,3, ...p) (A22 c) 


The stability criterion is the determinant of the coefficients of the Lagrangian multi pllerB in 
equations (A22) . 

numerical results .— numerical results for the critlcal-shear-e tress coefficient obtained in the 
present case are as follows : 


Approximation 

Lower limit 

Upper limit 

Firstj p*q«l 

12.50 

13.56 

Second ; p = q = 2 

13.06 

13.22 
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Thus, the average of the upper— limit and lower— limit results of 
the second approximation must be within 0.6 percent of the true 
buckling stress coefficient. 

Lower— limit results for the other types of buckling configurations 
were sufficient to indicate that the buckling configuration Just 
analyzed corresponds to the lowest buckling stress. Analysis of the 
second configuration (equation (Alb)) provided the lower— limit 
result k s = Ik . 29 with a sixth— order determinant; calculations for 

the third buckling pattern (equation (Ale)) gave the lower-limit 
value k 8 = lk.l with a fourth-order determinant. Since both of 

these values are higher than the upper— limit result for the first- 
buckling- configuration (equation (Ala)), it is evident that this 
first configuration is the governing one for the case of square bays. 

Although the same general method of analysis outlined in this 
appendix may be applied to the problem of the buckling of a continuous 
plate where the bays are not square, care must be taken to investigate 
all the possible buckling configurations; it may no longer be true 
that each bay will exhibit the same configuration as every other bay. 
Also, as the depth of the bays increases, the effect of continuity 
decreases, until, when b/a is infinite (see fig. 1), adjacent bays 
act as simply supported strips. 
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APPENDIX B 

DOUBLY INFINITE ARRAY OF SQUARE PANEIS 


The shear "buckling of the plate shown in figure 2 will be analyzed 
by a method closely analogous in general principle and specific detail 
to that of appendix A. Two possible buckling configurations will be 
investigated; as before, the governing configuration is that which 
provides the lowest buckling load. 

Boundary and continuity conditions, and deflection functions .— In 
this problem it is quite evident that the buckling configuration will 
be symmetrical about the midpoint of each bay; it is only necessary, 
therefore, to investigate the question of whether the direction of 
buckling Is the same or alternating from bay to bay. The two possible 
configurations may be represented by the following series, which satisfy 
the required continuity conditions term by term: 





u— b-H 


_ 

. + s 

+ s 

+ s 

\ 

b 

* 


+ s 

+ s 

+ s 



+ s 

-X 

+ s 

+ s 








\ \ , 2mrtx , 2nny \ \ , 2mjtx 2nny 

w = 2L 2_ ^ sin “b“ sin “TT + 2_ 2_ ^ 008 ~ b 008 “b 

m=l n=l m=0 n=0 


|S(0,J) - |5fb,y> 

ox ox 


Y (Bla) 


~( x ,o) = b) 
oy oy 
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b-b-H 


y, 

s 

+ s 

- s 

♦ 

b 

i 


+ s 

- s 

+ s 



- s 

\ 

+ s 

-s 








v 




00 

E 

n.3=l, 3i5, • • • 


asm 008 


mjtx 

b 


COB 


nity 

b 


00 


* E 



n=l,3,5. 


V 8in Tp Bin 


nrcy 

b 


> 


(Bib) 


|[(o,y) =-^B jy ) 

|=(x, 0 ) --| 2 (x,b) 


The boundary conditions of zero deflection at the supports must be 
introduced by meanB of Lagrangian multipliers. 

Stability crlterlons . - The first of the aforementioned buckling 
configurations (equations (Bla)) was the governing one; the remainder 
- of this analysis will be performed for that case. 

It may be seen that the deflection function chosen to represent the 
buckling configuration in the present problem is identical with that 
employed in the detailed analysis of appendix A. As a result of the 
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close si milar ity of the two problems, much of the analysis of appendix A 
may he applied to the present problem with only simple revisions. The 
constraining relationships are the same for the two cases except that in 
the present problem no condition of zero slope at the boundaries exists; 
therefore, the condition of equation (A 7 ) is omitted, and, in 
equations (A 9 ) and (A21), the multiplier is set equal to zero. 

A further simplification occurs in the present problem because of 
the symmetry about both diagonals of the buckling pattern within each 
bay. As a result it can be proved that, among the terms of deflection 
function w. 


a mn = a nm 
&mn = ^raa 


Consequently, in the application of equations (A9) and (A21) to the 
present case. 


a * 0 


and 


Xj = for i = j 


Hence the stability equations may be written directly from 
equations (Al 8 b) and (A22a) by employing the following relations: 


TJi = 0 
a = 0 

Xj = Hi for i = J 


First, for the lower— limit result, from equation (AlSb), the 
stability equations are 


00 



m=0 


if 

z 


+ y 

m=l 


m 


= 0 


(J = 1,2,3,... p) (B2) 
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For the upper^-limlt result, from equation (A23a), the stability 
equations become 

P p 

= 0 (J = 1,2,3,... p) (B3) 

m=0 m=l 

Numerical results .— The computed results for k a in the present 
problem are tabulated belov: 


Approximation 

Lover limit ' 

Upper limit 

First; p = 2 

11.07 

11.18 

Second; p = 3 

11.10 

11.10 


The lover— limit result for k 0 obtained from the other buckling 
configuration (equation (Bib)) is k a = 13.10; this result indicates 
that the first pattern is the governing one . 

Extension of this analysis to bays other than square is subject 
to the same qualification concerning the introduction of other possible 
buckling configurations as that stated in appendix A. 
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